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The metal-insulator transition in VO2 was investigated using the three-band Hubbard model, 
in which the degeneracy of the 3d orbitals, the on-site Coulomb and exchange interactions, and 
the effects of lattice distortion were considered. A new scenario on the phase transition is 
proposed, where the increase in energy level separation among the tig orbitals caused by the 
lattice distortion triggers an abrupt change in the electronic configuration in doubly occupied 
sites from an S = 1 Hund's coupling state to a spin S — state with much larger energy, 
and this strongly suppresses the charge fluctuation. Although the material is expected to be a 
Mott-Hubbard insulator in the insulating phase, the metal-to-insulator transition is not caused 
by an increase in relative strength of the Coulomb interaction against the electron hopping 
as in the usual Mott transition, but by the level splitting among the tig orbitals against the 
on-site exchange interaction. The metal-insulator transition in Ti2 03 can also be explained by 
the same scenario. Such a large change in the 3d orbital occupation at the phase transition can 
be detected by linear dichroic V 2p x-ray absorption measurements. 

KEYWORDS: VO2, metal-insulator transition, Mott-Hubbard insulator, exchange interaction, lattice dis- 
tortion, Ti 2 3 , V 2p XAS 



Vanadium dioxide undergoes a metal-insulator tran- 
sition (MIT) at T = 340 K. 1 ) Simultaneously, a lattice 
distortion from the high-temperature rutile (R) struc- 
ture takes place, and in the low-temperature monoclinic 
(Mi) phase, V ions pair and the unit cell doubles. Each 
V ion (nominally V 4+ with 3d 1 configuration) in the VO2 
crystal is surrounded by an octahedron consisting of six 
nearest-neighbor oxygen ions (see Fig. 1) and an octa- 
hedral crystal field splits the 3d levels into twofold e g 
and threefold t^g levels. The energies of the low-lying t2 g 
orbitals are further differentiated by a small Dih crys- 
tal field and the orbitals can be described as \a) = \xy), 

kl) = 7f{l^> + \ zx )} and K 2 > = Tfil^) " \ zx )} us- 
ing x-, y- and ^-coordinates locally defined for each V 
ion, as shown in Fig. 1. In the rutile phase, V ions form 
a linear chain along the c-axis, and oxygen octahedrons 
surrounding the V ions share their edges in each chain. 
Two kinds of chains having the shared edges along [110] 
and [110] directions are present. While the a orbitals 
have lobes along the chain and are strongly hybridized 
with their neighbors on the chain, the coupling among 
the 7T orbitals on the chain is weak. In the Mi phase, V 
ions on each chain form pairs and the V-V bonds of the 
pairs tilt from the c-axis. 

A number of theories based on Peierls, 2-4 ' Spin- 
Peierls 5 ) and Mott-Hubbard 6 ' 7 -* mechanisms have been 
proposed for the MIT. Since V ions form pairs in the 
Mi phase, the Peierls mechanism appears to be a good 
candidate as an explanation for the MIT. However, ex- 
perimental results obtained from studies on the Cr-doped 
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Fig. 1. Crystal structure of V02- The large and small circles de- 
note oxygen and V ions, respectively. The unit cell in the mon- 
oclinic (Mi) phase is indicated by solid lines and those in the 
rutile (R) phase by dashed lines. The arrows indicate exagger- 
ated displacement of the V ions in the Mi phase. 



alloy Vi- x Ct x 02 S ' > and on VO2 under uniaxial stress 9 - 1 
are inconsistent with such an interpretation. In addition 
to the Mi and R phases, insulating monoclinic (M2) and 
triclinic (T) phases at the temperature between the Mi 
and R phases are found in Cr-doped samples and VO2 
under pressure. In the M2 phase, V ions on one-half of 
the chains pair without the tilting of their V-V bonds, 
whereas those on the other half form zigzag chains with- 
out pairing. These zigzag chains with equally spaced V 
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ions behave magnetically as S = 1/2 Hciscnberg chains 
with J ~ 300 K and are considered to be Mott-Hubbard 
insulators. Furthermore, a continuous M 2 -to-Mi transi- 
tion through the intermediate T phase is observed as 
temperature decreases. Since these two insulating phases 
appear by doping a small amount of Cr or uniaxial pres- 
sure, all the Mi , M 2 and T phases are classified as Mott- 
Hubbard and not band insulators. 10 ' 

The purpose of this letter is to discuss the mechanism 
of the MIT in VO2 using the three-band Hubbard model 
on the basis of the many-body theory, particularly in re- 
lation to the on-site exchange interaction and the level 
splitting among the t 2g orbitals caused by the lattice dis- 
tortion. 

As a model for the V ion chain along the rutile c-axis 
in VO2, a linear chain cluster consisting of six V ions 
with the periodic condition is assumed and the three t 2g 
orbitals a, 71T and 7r2 are considered in each V ion. The 
Hamiltonian assumed in the model is written as 



exchange interactions: 



H, 



clcc 



H cry + #hyb + H\ 



(i) 



In the band structure calculations, the energy band aris- 
ing from the n orbitals is shifted to higher energy due 
to the pairing and zigzag distortions of the V ion chains 
in the insulating phases. 2,4 ' To consider such effects of 
the lattice distortions in the Mi and M 2 phases, a low- 
symmetry crystal field iJ C ry which splits the a and n 
orbital levels is introduced and is described as 
D 



H c 



■ ^2 (n, j7r i + n li7r2 - 2m >(T ) . 



(2) 



where D is the energy difference between the a and n 
orbitals and n ijjrl is the number operator for an orbital 
m(= a, 7rl or n2) at the i-ih V ion site. 

-ffhyb represents electron hopping between the nearest 
neighbor V ion sites and is written as 

-ffhyb = ^2 V™ ( a i,m,s fl i+l,m,s + a i+l,m,s a i,m,s) i (3) 

where a\ m s denotes the creation operator for an elec- 
tron on an orbital m with a spin quantum number s at 
the i-th V ion site and V m is the hopping integral. For 
simplicity, the variations in the values of the hopping in- 
tegrals arising from the pairing and zigzag distortions of 
the chains in the insulating phases are not considered. 
However, even if these effects are included, there is no 
qualitative change in the results obtained. Thus, D is the 
only essential parameter of the distortion for the mech- 
anism of the MIT, and the Peierls and Mott-Hubbard 
mechanisms triggered by the change in the strength of 
the hybridization cannot be the main cause of the MIT 
within the present model. 

The last term H mt stands for the on-site Coulomb and 
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where U and U' are the Coulomb repulsion energies be- 
tween electrons on the same and different orbitals, re- 
spectively, and J is the exchange interaction energy; 
there exists a relation U = U' + 2J among them. The 
values for U, U' and J are determined from high-energy 
spectroscopic data. Note that although only the V 3d or- 
bitals are included in the present model, the hybridiza- 
tion between the 3d and neighboring oxygen 2p orbitals 
is strong in the actual system. To consider the effects of 
the 3d-2p mixing, the parameter values were adjusted 
to reproduce a low-energy level diagram obtained from 
a cluster model, in which both the V 3c? and O 2p 
orbitals are included. 11 ' For the hopping integrals be- 
tween the t 2g orbitals, not only the direct d-d hopping 
but also those via neighboring oxygen 2p orbitals is also 
considered within the second-order perturbation theory: 
V a = (3dda+dd8)/{AN^), V nl = {ddir-pdir 2 / 'A) /N n and 
Vtt2 = (ddS + pdn 2 / A) / 'N„ . Here, dda, ddn and ddS are 
the Slater-Koster parameters for the 3d orbitals between 
the nearest neighbor V ions, pdir is that between the 3d 
and 2p orbitals, A is the charge transfer energy between 
the 3d and 2p orbitals and iV^ = 1 + 4(pc?7r/A) 2 . For 
the values of A and pdir, those obtained from analysis of 
the high-energy spectroscopy experiments are used, and 
for dda, ddn and ddS, Harrison's values are adopted. 12 ' 
The parameters adopted in the model are U' — 1.33 eV, 
J = 0.65 eV, V a = -0.25 eV, V vl = -0.04 eV and 
K-2 = 0.22 eV. The ground-state and low- lying excited- 
states eigenvalues and eigenvectors were calculated by 
numerical diagonalization of H c i cc using the Lanczos 
method. 

A possible scenario on the MIT is as follows. When the 
value of D is small, the electronic configuration an (one 
of the two electrons occupies the a orbital and the other 
the 7r2 orbital) with total spin S = 1 for the doubly occu- 
pied site is favorable because of Hund's first rule, and this 
causes nearest-neighbor ferromagnetical antifcrro-orbital 
correlation as a]— 7r|. On the other hand, for a large D 
compared to J, the S = aa configuration is favorable 
for the doubly occupied site and this induces nearest- 
neighbor antifcrromagnctical ferro-orbital correlation as 
a]-a[. Because on-site electron repulsion energy for the 
S = 1 an configuration is V — J = 0.7 eV whereas that 
for the S = aa configuration is U = 2.6 eV, charge flue- 
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Fig. 2. D dependence of the ir2 and a orbital occupation per V 
ion (a) and weight of the d 1 d 1 d 1 d 1 d 1 d 1 configuration (b) in the 
ground state of the V ion chain model. For the weight, results 
obtained with various values of J are depicted. 



tuations c? 1 ;^ 1 — > d°;d 2 should be strongly suppressed 
for large D. 

Figure 2 shows D dependence of the n2 and a or- 
bital occupation per V ion (a) and the weight of the 
d 1 d 1 d 1 d 1 d 1 d 1 configuration (b) in the ground state of 
the V ion chain cluster model. To determine the effects 
of J on the charge fluctuation, in Fig. 2(b), the weights 
obtained with various values of J = 0, 0.32, 0.65 and 
0.97 eV are shown. As anticipated in the above discus- 
sion, for D < 0.1 eV, the ratio of the orbital occupation 
is no- : n w = 1 : 1, which implies the presence of the 
S = 1 an doubly occupied sites and erf— 7r| correlation 
in the ground state, and for D > 0.34 eV, the ratio is 
= 1:0, which indicates the presence of the 
S = aa doubly occupied sites and ctT _(J I correlation. 
Since we find other states with intermediate orbital oc- 
cupations n a : n x = 2 : 1 and n a : — 5 : 1 in 
0.10 < D < 0.34, the variation from the state with the 
5* = 1 an to S = aa doubly occupied sites with in- 
creasing D is probably continuous in the infinite-sites 
limit. In Fig. 2(b), we can find an approximately 60% 
increase in the weight of the d 1 d 1 d 1 d 1 d 1 d 1 configuration 
from the ground state around D = eV to that around 
D = 0.4 eV with J = 0.65 eV. This indicates a large 
change in the charge fluctuation between the two states. 
As discussed earlier, this change in the charge fluctuation 
originates from the difference in the electron repulsion 
energy of the doubly occupied sites: U' — J for the S = 1 
an configuration and U(= U' + 2J) for the S = aa 
configuration. Indeed, both the difference in the weights 
and the values of D where the change in the configura- 
tions take place increase with increasing J, showing that 
the on-site exchange interaction is essential for the MIT. 

The two kinds of electronic configurations discussed 
above are analogous to the so-called high-spin and low- 
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Fig. 3. Optical conductivity spectra calculated with D = eV 
(a) and D = 0.45 eV (b). 



spin states in a transition metal ion embedded in an oc- 
tahedral crystal field scaled by the energy difference be- 
tween the e g and ti g orbitals lODq. If the value of lODq 
is small with respect to the intra-atomic exchange inter- 
action energy, the 3c? electronic state still has the same 
total spin as expected from Hund's first rule and is called 
a high-spin state, whereas for the crystal field superior 
to the exchange interaction, an electronic state with all 
electrons in low-lying t2 g orbitals having smaller total 
spin is favorable and is called a low-spin state. From this 
viewpoint, the phase transition can be interpreted as a 
"high-spin" 5* = 1 an state to "low-spin" S — aa state 
transition in doubly occupied sites caused by a change in 
the crystal field D. 

Optical conductivity spectra were calculated to deter- 
mine the difference in electronic transport properties be- 
tween the ground states with the low-spin and high-spin 
doubly occupied sites corresponding to the insulating 
and metallic phases in real material, respectively. The 
method described in rcf. 13 is employed for the calcula- 
tion with the six V ion cluster model. Figure 3 shows the 
optical conductivity spectra as a function of the photon 
energy uu calculated with D = eV (a) and D = 0.45 eV 
(b) . The spectra were convoluted with the Gaussian dis- 
tribution curve with two different broadening widths 0.02 
eV and 0.4 eV. As seen in Fig. 3(a), a large portion of the 
intensity ^95% is concentrated at the Drude peak, which 
is positioned at u — 0, and this clearly shows the metallic 
feature of the ground state with the S = 1 an configu- 
ration in the doubly occupied sites. On the other hand, 
the spectrum for the ground state with the S = aa 
configuration in Fig. 3(b) has a massive peak at ~2 eV 
and the intensity of the Drude peak is reduced to only 
~ 5% from that in Fig. 3(a). Although the weak Drude 
peak still remains, probably due to the finite size effect, 
such a large reduction in its intensity indicates the local- 
ization of electrons in the ground state with the low-spin 
doubly occupied sites. 

We have discussed the variation in the electronic struc- 
ture caused by the monoclinic lattice distortion scaled by 
D. In the real material, the first-order MIT takes place 
as temperature decreases accompanied with the lattice 
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Fig. 4. Temperature dependence of D which minimizes the free 
energy F(T; D) (a) and the expectation value of n CT (b) obtained 
with K = 1.17 (eV)- 1 and various values of Dq. 

distortion. To demonstrate how the change in the elec- 
tronic structure as a function of D can be related to the 
temperature- induced MIT in VO2, let us consider a sim- 
ple model for the lattice elastic energy, where the energy 
per V ion site is a quadric function of D: -§-(-D — -Do) 2 - 
The free energy F(T; D) of the system then can be writ- 
ten as 

F(T;D) = -k B T\ogZ clcc + N sitc ^(D-D ) 2 , 

where Z e \ ec is the partition function of H e \ ec defined in 
cq. (1) and the second term is the lattice elastic energy 
and -/Vgitc = 6 is the number of V ion sites in the model. 
Figure 4 shows values of D which minimize the free en- 
ergy F{T;D) for each T calculated with K = 1.17 (eV)" 1 
and D = -0.10, -0.11, -0.12, -0.13 and -0.14 (eV). In 
Fig. 4(b), the expectation values of the a orbital occupa- 
tion are also depicted. In Fig. 4(a), a clear discontinuous 
change in D as a function of T can be seen. With this 
jump in D, the electronic state also changes from the lo- 
calized state with the S — era sites to the delocalized 
one with the S — 1 a ir sites, which is evident from Fig. 
4(b). This temperature-induced abrupt transition, both 
in the lattice distortion and the electronic state, which 
is consistent with the experimental MIT, can be under- 
stood as follows. For the state with S = 1 an doubly 
occupied sites, the excitation originated from its spin, 
and orbital degrees of freedom positioned densely within 
the low-energy region < 0.1 eV lower the free energy at 
high temperatures. On the other hand, such an energy 
lowering by the entropy term is minor for the state with 
the S = aa sites, since the state has only spin excita- 
tions. 

If there is a large change in the electronic configura- 
tion between the metallic and insulating phases, as is 
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Fig. 5. V 2p XAS spectra for two different directions of the polar- 
ization vectors e, which are parallel (solid line) and perpendicular 
(dashed line) to the rutile c-axis, calculated with D = eV (a) 
and D = 0.35 eV (b). The difference spectra are shown at the 
bottom in each panel. 

predicted in the present theory, such a difference can 
be detected by measuring the linear dichroism in the 
V 2p x-ray absorption spectrum (XAS), since the spec- 
trum is very sensitive to the 3d orbital occupation. To 
show this, the spectra were calculated using a cluster 
model with four V ions, where, in addition to the ti g 
orbitals, 11,14 ) the e g orbitals arc considered for each V 
ion site. While the ground state properties are almost 
unchanged, inclusion of the e g orbitals is important for 
the final state of the spectra. In the final state, the 2p 
core hole of photoexcited site, and the Coulomb and ex- 
change interactions between the 3d electrons and the 2p 
core- hole are further taken into account. Figure 5 shows 
the theoretical spectra obtained with D = 0.0 eV (a) and 
D = 0.35 eV (b) corresponding to those for the metal- 
lic and insulating phases, respectively. The spectra are 
depicted for two different directions of the polarization 
vector e of the light: parallel and perpendicular to the 
rutile c-axis. The difference spectra are also shown at 
the bottom in each panel. Because of the selection rule 
of the 2p — > 3d dipole transition, the transition probabil- 
ity to the a orbital is larger than that to the tt orbitals 
for e//c r , whereas that to the 7r orbitals is larger for the 
£ _L c r . Indeed, the integrated intensity of the difference 
spectra for the "metallic" ground state with the orbital 
occupation n a : n n =l : 1 and the "insulating" ground 
state with n a : n n =l : have mutually opposite signs, 
in addition to the large differences in the spectral shape 
itself. 

An MIT in Ti2C>3 can also be explained by the same 
mechanism proposed here. Ti2C>3 exhibits a gradual MIT 
around 400-600 K and has the corundum structure, 
where Ti ions form pairs along the c/j-axis. 1 ' Because 
of a small trigonal field, the ti g levels split into low- 
lying ai g and twofold levels. At high temperature, 
the electronic configuration of the doubly occupied sites 
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is expected to be S — 1 ai g e g 7T . Although the symmetry 
of the crystal is not changed, the Ti-Ti bonds of Ti ion 
pairs become shorter and the energy separation of the e g * 
and ai g orbitals increases as temperature decreases, and 
this causes the high-spin to low-spin (5 = ai g ai g con- 
figuration) transition at the doubly occupied sites. Since 
the ai s and e g n orbitals are hybridized through inter- 
pair hopping integrals, the transition is not as steep as 
we have seen in the V0 2 calculations. However, a large 
change in the orbital occupation from n alg : n eg ~ 1 : 1 
to n a i g : n eg ~ 1 : is expected. 

In conclusion, the MIT in VO2 has been discussed us- 
ing the three-band Hubbard model in connection with 
the lattice distortion on the basis of the many-body the- 
ory. A novel mechanism on the MIT is proposed. While 
the material is considered to be a Mott-Hubbard insu- 
lator in the insulating phases, the MIT is not induced 
by an increase in the relative strength of the Coulomb 
interaction against the electron hopping, as in the usual 
Mott transition, but by the level splitting between the a 
and 7r orbitals against the on-site exchange interaction, 
where the charge fluctuation is suppressed by the change 
in electron configuration in doubly occupied sites from a 
high-spin (S = 1) to a low-spin (S = 0) state with much 



larger interaction energy. 
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